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Abstract 

Starting from the T-Q equations of the open spin-1 XXZ quantum spin chain with 
general integrable boundary terms, for values of the boundary parameters which sat- 
isfy a certain constraint, we derive a set of nonlinear integral equations (NLIEs) for 
the inhomogeneous open spin-1 XXZ chain. By taking the continuum limit of these 
NLIEs, and working in analogy with the open spin-1 XXZ chain with diagonal bound- 
ary terms, we compute the boundary and the Casimir energies of the corresponding 
super symmetric sine-Gordon (SSG) model. We also present an analytical result for the 
effective central charge in the ultraviolet (UV) limit. 
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1 Introduction 



Due to applications in statistical mechanics and condensed matter physics, one dimensional 
spin systems have attracted much interests. These models have been subjected to intensive 
studies over the years. One such model in particular is the XXZ quantum spin chain with 
boundaries. For example, various aspects of the open spin- 1/2 XXZ chain have been actively 
investigated [1]-[15]. Among others, its finite size properties of ground and excited states 
have received much attention over the years. This lattice model has also been used to 
study the ground and excited states of the corresponding field theoretical model i.e., the 
sine-Gordon (SG) model with boundaries [16], let it be with Dirichlet or general integrable 
boundary conditions [14], [17]- [28]. As for the periodic case [29]-[31], NLIEs have also been 
derived from the corresponding Bethe ansatz solutions of the open spin- 1/2 XXZ chains with 
nondiagonal boundary terms [9]- [13], [15] and used in such investigations [14, 22, 26, 27, 28], 
in particular, in studies involving the crucial ultraviolet (UV) and infrared (IR) limits. 

The spin-1 XXZ chain which is related to the supersymmetric sine-Gordon (SSG) model 
[32] -[46] has also gone through similar studies over the years. The NLIEs have been obtained 
for cases with the periodic boundary conditions [45] and Dirichlet boundary conditions [46]. 
They have been used to compute quantities such as S matrices (both bulk and boundary), 
bulk and boundary energies, central charges and conformal dimensions. However, as pointed 
out in [46], since the ground state of the critical spin-1 XXZ chain is described by a sea of 
approximate "two-strings", the usual method [29]- [31] of deriving these NLIEs from Bethe 
ansatz equations and the counting functions does not seem to work in this case. Fortunately, 
NLIEs can be derived [47, 48] rather by exploiting the analyticity of the transfer matrix 
eigenvalues described by the T-Q equations of the model. Such a method have been used to 
derive the desired NLIEs describing SSG models with periodic boundary conditions [45] and 
the Dirichlet boundary conditions [46] and used to calculate finite-size properties of ground 
and excited states. 

Our motivation for carrying out this work comes from these works, where spin-1 XXZ 
chains have been associated with SSG models. Moreover, solutions for the open spin-1 XXZ 
chain with general integrable boundary terms [49] have been recently proposed [50]. Hence, 
using this solution, one will be able to derive a set of NLIEs for the inhomogeneous open 
spin-1 XXZ chain which at the continuum limit should describe a SSG model. And since 
such NLIEs from open spin-1 XXZ chain with general integrable boundary terms have not 
been derived, we plan to address this issue in this paper. We shall only consider the ground 
state in this paper. 

The outline of the article is as follows. In section 2, the Hamiltonian and T-Q equations 
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for the open spin-1 XXZ spin chain with general integrable boundary terms are reviewed. 
This is followed by the derivation of a set of NLIEs for the inhomogeneous open spin-1 XXZ 
chain in section 3. We then give a brief review of the SSG model following [41]. Proceeding 
by analogy with the SG model with general integrable boundary conditions [14, 27] and the 
SSG model with Dirichlet boundary conditions [46], we consider the continuum limit of the 
NLIEs in section 4. In section 5, we compute the boundary and Casimir energies for the 
ground state. We analyze the UV limit of the Casimir energy and give an analytical result 
for the effective central charge in section 6. Finally, we conclude the paper with a brief 
discussion of our results and some open problems in section 7. 

2 The open spin-1 XXZ chain 

2.1 Hamiltonian of the open spin-1 XXZ chain 

In this section, we first begin by reviewing the Hamiltonian of the open spin-1 XXZ chain. 
We adopt the notations used in [50]. 

JV-1 

U=Y, H n,n + l + H b . (2.1) 

71=1 

where i? n ,n+i represents the bulk terms. Explicitly, these terms are given by [51], 

Hn,n + 1 = On' + 2 sh 2 T) [< + {S^f + (^ +1 ) 2 - «) 2 ] 

- 4sh 2 (|)(^< + <^) , (2.2) 

where 

&n = S n ■ S n+1 , = S^S^ +1 + S^S y n+1 , a z n = S^S'^ +1 , (2.3) 

and S are the su(2) spin-1 generators. represents the boundary terms which have the 
following form (see e.g., [49]) 

H b = ai {Slf + a 2 Sl + a 3 (S+) 2 + a 4 (S^) 2 + a 5 S+ S* + a & S{ 

+ a 7 St Sf + a 8 Si S z x + (a, b j and loJV), (2.4) 

where S ± = S x ±iS y . The coefficients {aj} of the boundary terms at site 1 are functions of 
the boundary parameters (a_, (3_,9_) and the bulk anisotropy parameter 77. They are given 
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by, 



where 



a 2 
a 3 

0,4 

a 5 
a 6 
a 7 
a 8 



-a (ch 2a_ — ch 2/3_ + ch 77) sh 2rj sh 77 , 
-a sh 2a_ sh 2/3_ sh 2i] , 



— -a e 26 sh2r]sh.r] : 
8 

— -a e~ 2e - sh.2rjsh.ri, 
8 



a o e ~ yzh /3_ sh a_ ch - + ch a_ sh /3_ sh - J sh 77 ch 2 77 , 

aoe~ e_ ( ch /3_ sh a_ ch — h ch a_ sh /3_ sh - j sh 77 ch 2 77 , 
V 2 2 / 
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-aoe " ( ch (3_ sh a_ ch — — ch a_ sh /3_ sh — ) sh 77 ch 2 77 . 



-a e" 



( 



77 \ 3 

ch /3_ sh a_ ch '4 — ch a_ sh /3_ sh - J sh 77 ch 2 77 , 



77 
2 



a 



sh(a_ - 1) sh(a„ + i) ch(/3_ - i) ch(/3_ + ±) } 



-1 



(2.5) 



(2.6) 



Similarly, the coefficients {hi} of the boundary terms at site N which are functions of the 
boundary parameters (a + , 9+) and 77, are given by the following correspondence, 



(2.7) 



2.2 The T-Q equations for the inhomogeneous open spin-1 XXZ 
chain 

Next, we review the T-Q equations for the open spin-1 XXZ chain with general integrable 
boundary terms given in [50]. However, we note that the solution holds only when the 
boundary parameters (a± , (3± , 9±j obey the following constraint [11]- [13], [50], 

a_ + /3_ + a + + /3 + = ±(9_-9 + )+r]k, (2.8) 

where k is an odd integer. A convenient redefinition of bulk and boundary parameters can 
be adopted [14]: 

77 = i/i , a± — i/ia± , j3± = fib± , 9± = i[j,c± , (2.9) 
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where /i , a± , b± , c± are all real, with < ji < |*. With the above redefinitions, the constraint 
relation (2.8) assumes the following pair of real contraints: 

a_ + a + = ±|c_ — c + \ + k , 

b_ + b + = 0. (2.10) 

In this paper, we shall consider only even N case. Two relevant commuting transfer matrices 
for spin-1 XXZ chain are Ti(u) with a spin-^ (two-dimensional) auxiliary space, and T 2 {u) 
with a spin-1 (three-dimensional) auxiliary space. The corresponding eigenvalues obey T — Q 
equations found in [50]: A\(u) which represents the eigenvalues of Ti(u) can be written as 
(following [46] and adopting the notations used there) 

Ai(«) = h(u) + l 2 (u), 

Uu) = sinh(2 M + tfi)B i+ \u)<j){u + in) ^ , 

Q{u) 

Uu) = sinh(2^ - in)B^\u)4>(u - ^) ^'" + , , (2.11) 

Q(u) 

where 

<j>{u) = sinh 7V (M-A)sinh JV (M + A), 

~( + \, i/j,A + i/j,A^ i/j,B + i/j,B_ 

B y >(u) = smh(u ± — - — ) sinn(M ± — - — ) cosh(-u =F — ^ — ) cosh(-u =F — ^ — ), 

Zi Zj Zi Zi 

M 

Q(u) = |^sinh(-u — v k ) sinh(-u + v k ) . (2-12) 

k=l 

We have redefined the lattice boundary parameters as A± = 2a± — 1 , B± = 2ib± + l.In 
addition, A is the inhomogeneity parameter, N is the number of spins and M = N — | + | 
represents the number of Bethe roots which are also the zeros v k of Q(u). In this paper, we 
shall consider the case k — 1. For this particular model of the XXZ chain, one generally 
needs to consider two groups of transfer matrix eigenvalues, labelled as A^' 1 )( ± )(-u) in [50], 
to obtain all 3^ energy eigenvalues. Since we are interested only in the ground state, we 
restrict our analysis to only one of them that contains the ground state, namely A^' 1 )( _ )(-u). 
Readers are urged to refer to [50] for details on this. Next, using the fusion relation, one can 
write the eigenvalues of T 2 (u), A 2 (u) as (see e.g. [50]) 

A 2 («) =A 1 (u- |) A 1 (u + |) - f(u) (2.13) 



1 Refer to section 4 of [46] on IR limit for details on this bound. 
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where 

f{u) = <f ) (u+^) ( j)(u-^)smh(2u-2i^smh(2u + 2i^B^\u + ^)B^\u- i -^) 
2 2 2 2 

= h(u+ l p 2 (u-^) (2.14) 
Using (2.14), (2.13) can thus be written as 

Hu) = / 2 («-^)Z2(«+^) + Zi(«-^)Z 2 («+^) + /i(«-^)Zi(«+^) 

= sinh(2w)A 2 (w) (2.15) 

where 

A 2 (it) = smh(2u-2t^-\u-^)B^\u + ^) ( f>(u-^) ( f>(u 



2 /ry 2 2' Q(u-f) 

+ smh(2.)S( )(«- y )B( )(u + T )^(u- T)< Ku + T ) q(u _ ¥) 

+ sinh(2. + 2 ¥)5 W(n + - | + ^)0(. + f 

= Ai(m) + A 2 (m) + A 3 (m) (2.16) 
From (2.15) and (2.16), we also have 

-mW — • i/ \ 5 A 2 (MJ — . , . , <A 3 (WJ 



sinh(2u) ~ sinh(2-u) ' sinh(2-u) 

(2.17) 

One can now define the auxiliary functions b(u) and b(u) by 

b(u) = hH±^± , b(u) = b(-u) = hH±^ (2.18) 
X 3 (u) Ai(u) 

We note that 6(it) is the complex conjugate of b(u) for real u. Using (2.11), (2.17) and (2.18), 
we obtain 



sinh(2w + 2in) <j>( u + f)<j>(u + ?f) B(+)(u - f + §) Q(u - ^) 



We also note that Ai(u) does not have zeros near the real axis, except for a simple zero at 
the origin. Following [46], one can remove this root by defining, 

A lM ^ „ 20 ) 



where k(u) is any function whose only real root is a simple zero at the origin, that is 
k(0) = , k'(0) 7^ 0. The prime denotes differentiation with respect to u. In terms of Ai(u), 
(2.19) can be compactly written as, 

b(u) = C b (u) K^u - it) Q ^ + f^ (2.21) 



k{u - f )<j>{u - | ) gR^+ji) 

sinh(2« + 2i^)0(« + f )0(m + f*) - f)B(+)(u + f ) 



where 

w = -^/?:^ 2 rr^:, , . ^ ( 2 - 22 ) 

Next, defining 

£(u) = l + &(«), £(u) = l + 6(u), (2.23) 

and using (2.16) and (2.18), we have the following expressions for A 2 (u), 

A 2 (u) = A 3 («)S(u) 

= Ai(w)S(u) (2.24) 

As for Ai(u) (see (2.20)), A 2 (w) also has a root at the origin which can be removed by 
defining 

*„>-^>. „ 25 ) 

Using (2.11), (2.17) and (2.24), (2.25) becomes 

C 

Q(u-f) 
in 

Q(u+i) 



Mu) = t + (u) Q ^ +3 lh (u) 



= M*) t£ ■ l} B(u) (2.26) 



where 

= ■mi.^i p,, _ | )J3 , T)(11 + i )0( „ T ^ T | ( . (2 . 27) 

z III 

We next define the last two auxiliary functions y(u) and Y(u) as follows, 

. , sinh(2-u)A 2 (w) . . , , , , 

y{u) = y f ^ 21 ; , y(u) = l + y(u), (2.28) 
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2.3 Bethe roots and parameter regions for ground state 



In this paper, we are primarily interested in studying the finite size properties of the ground 
state. Hence, the distribution of Bethe roots, {vk} for ground state of the open spin-1 
XXZ chain described by the Hamiltonian (2.1), with the boundary parameters satisfying the 
constraint (2.8), must first be given. The ground state is described by sea of approximate 
"two-strings", Vk = %k ± Wk (where Xk and y& are real), the schematic depiction of which 
is given in figure 1 below. One can verify numerically that the imaginary parts yk satisfy 



< y k - £ < i. 






in/2 




-in/1 



Figure 1: Zeros of Q(u) 

From numerical investigations, the regions in parameter space A± which yield the above 
form of Bethe roots for the ground state can be divided in the following way. 2 



I 
II 
III 
IV 



1 < A ± < T 
KA + <f & 

1 - f < A ± < - 

l-^<A + <- 



ft 



& 1< A_ < ^ 



(2.29) 



In addition, due to (2.10), the parameters B± satisfy 

B+ + B_ = 2 



(2.30) 



In subsequent sections, following [46], we will work instead with a shifted function, defined 
by q(u) = Q(u + y). This is due to the fact that for ji — > 0, Q(u) can have zeros near 
the real axis (refer to figure 1). On the other hand, q(u) does not suffer from such features. 
Refer to the figure 2 that gives schematic picture of zeros of q(u). 





>(7T - m)/2 




-»(tt - m)/2 



Figure 2: Zeros of q(u) 



2 Such investigations have been given in [14] for the corresponding spin-1/2 XXZ chain. Readers are urged 
to refer to the reference for further details on such numerical investigations. 



7 



3 The NLIEs for the inhomogeneous open spin-1 XXZ 
chain 



Utilizing the analyticity of 1iiA 2 (m) near the real axis, we have the following from Cauchy's 
theorem, 

= j> du [\nA 2 (u)]"e iku (3.1) 

where the contour C is chosen as in figure 3 below, e is small and positive, such that the 
max {y k } - f < e. 



Ci 



it 



— l€ 



c 2 



Figure 3: Integration contour 



After using (2.26), (3.1) can be written as 
= 



du {lnt_(u)} " e lku + / du <J In 
Ci Jc\ 



iku 



+ I du {\nB(u)}"e iku + 



Ci 



/ du {lat + (u)}" i 
Jc 2 



iku 



+ / du { In 

ic 2 



/ du {\nB(u)}"e iku , (3.2) 
•/c 2 



As in [46], we define the following Fourier transforms along C 2 and C 1: 



Lf"(k)= du [\nf(u)]"e zku , Cf"(k) = du [In f (u)}" e lku , 
>c 2 Jc\ 



(3.3) 



respectively. Exploiting the periodicity (to make the imaginary part of the argument nega- 
tive), 

q{u) = q{u — m), u G C±, and q(u + iu) = q{u + ifj, — in), u G C 2 . (3.4) 
we obtain the following 



/ 

J a 
la 



du < In 



du < In 



q(u - 



3iu , ITT 

2 2 



r' A " = Z,f (A;) (Vt-i)* - e (!-f)* 

= L?(fc) (e^-i) fe - e (f-t) fe ) , (3.5) 



Aku 



Consequently, after using (3.3) for t±(u), B(u) and B(u) together with (3.5), (3.2) becomes 

,(*-¥)* _ e (f-f ) k _ e (¥ -f )* + e (i-f = LB"(k) + + Lt i l(k) + Ct i L(k), 

(3.6) 

Also from (2.21), we have the following for the Fourier transform of b(u), 



Lb"{k) = e~ 2 LA'{(k) + Lq"(k) 
In addition, from (2.13) and (2.20), we obtain 3 



e(¥-f>_ e (f- 



(e^r + e - ^) LA"(k) = LY"(k) + Lf"(k) - (e^ + e"^) £«"(£;) + e^27rA; . 



(3.7) 



(3.8) 



which can be used together with (3.6) to obtain the following NLIE in Fourier space for b{u), 
LV'{k) = -G(k) l3 l (k)+jCB''(k)j+G 2 (k)LY } '(k) + C(k) (3.9) 

where 

sinh ((it - 3/i)|) 



G(k) = 



G 2 (k) = 



2 cosh f sinh ((tt- 2//) |) ' 



fj,k 

e 2 



e 2 + e 2 



(3.10) 
(3.11) 



C(k) = -G{k){Lt\{k) + CfL(k)) + G 2 {k) Lf"(k) + LC£(k) - e-^LK"(k) 
2irk 



+ 



fik fik_ 

e 2 + e 2 



(3.12) 



where the expressions for CV(A;), C*i(A;), and C^ ± (k) (evaluated using (2. 14), (2. 22) and 
(2.27)) are given below, 



C N (k) = 2irN^(k)(e lAk + e- lAk )[-G(k)(e~^+e-'^-e- k ^-e 



) 



+ G 2 {k){e~— +e 
d(k) = 2m/j 2 (k)[-G(k)(e 



2 ) + (e 2 - e 



+ 



)] , (3.13) 
) + G 2 (k)(e- k ^ + e k ^)] + 2n(kG(k) - ^ 2 {k)e k ^ 

(3.14) 



C & Jk) = -27r^(A;)|s_e"("^ + ") fc + S+e "("^ + ") fe - S ^ — " - „■_< — 



1 + e fc ^ + e~ fc ^ 



(- 



e 2 + e 2 



G(ib)(e^ + e-^)) 



(3.15) 



3 We remark that the term 2irk in (3.8) is obtained by evaluating the integral § c du [In n(u)\" e lku . 



where s± = sgn(A ± ), ^(/c) = i_ k -*k and = — ^4ff and we have used the following 

1 e l-e ~2~ 

identities (see also [46]), 

/ — [lnsinh( M -^)fe lfeu = e" fc( ™VW, (3.16) 

where n is an integer such that < 3?e(a — nir) < ir, and 

/ — [lnsinh(2 M )]"e ifcu = ^2(A;). (3.17) 
Jc 2 27r 

Evaluation of C(k) from (3.12) thus yields the following, 

C(k) = 2^hl N ( ^ + 6 ~ lAk ) g + sirih ((^ - + sinh ((tt - 

U n \ \ 2 cosh f J 2cosh(f )sinh((7r-2/x)|) 

sinh + sinh (|gg-) cosh ^ sinh ((3/x - tt)|) \ 

2 cosh(f ) sinh ((tt - 2//)§) cosh f sinh ((2/x - vr)|) J ' 

The third NLIE involves 

Ly"(k) = Lk'^k) + [ du [lnsinh(2w)]V fcM + Ltf'{h) - Lf'(k) . (3.19) 
Jc 2 

To this end, we proceed with the evaluation of LA' 2 \k). From (3.1) and (3.2), and through 
the elimination of Lq"(k), we have 

(e^ + e -f j Lk'^ik) = e~^LB"(k) - e^CB"(k) + e^L^A;) - e^£?(A;) . (3.20) 

Solving for LA'^k) from (3.20), and using the result in (3.19), we obtain 

L^'(k) = -G(-k) CB"(k) + G 2 (k) LB"{k) + C y {k) , (3.21) 

where 

C y (k) = ^ Lt/ m-e^Lt'L(k) + 2 ^ (A;) + _ (3 22) 
e 2 + e 2 

which after some manipulation yields the following, 

C tf (A;) = AixkG 2 (-k) . (3.23) 

Thus, (3.9) and (3.21) represent the NLIEs of the inhomogeneous lattice model of the spin-1 
XXZ chain with general integrable boundary terms in Fourier space. For further computation 
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sinh — 
cosh^ 



in subsequent sections, we express these NLIEs in coordinate space. By integrating twice, 
we obtain 

/co pco 
du' G(u — u — ie) ln(l + b(u' + ie)) - j du' G(u - v! + ie) ln(l + b(u' - ie)) 
-CO J —CO 

/CO I c 

du' G 2 (u -u' + ie) ln(l + y(u' - ie)) + ^iVtan" 1 - 

+ i du' R(u') + Cm , 

/OO /'OO 
du' G(u - u' - ie) ln(l + b(u' + ie)) + / du' G(u - u' + ie) ln(l + b(u - ie)) 
■CO J —CO 

/CO 
du' G 2 {u' -u + ie) ln(l + y(u' + ie)) - i2A^tan _1 
-OO 

i \ du R{u) — Cm , 
Jo 



sinh — 
m_ 

cosh^ 



— i 



\ny(u) = I du' G 2 (u — u' + ie) ln(l + b(u' — ie)) + / ciw' G r 2( , w / — u + ie) ln(l + b(u' + ie)) 



+ Ani I du' G 2 (-u') + Cyin . (3.24) 
Jo 



where R(u) refers to the Fourier transform of R{k) which is given below, 

s + sinh ((7r - + s_ sinh ((it - 



= 2tt< 



2cosh(f )sinh((7r-2/i)|) 



sinh (|^g+) + sinh (IfiB.) cosh(f ) sinh ((3^ - tt)|) 
2cosh(f )sinh((7r-2/i)|) cosh(f )sinh((2/i-vr)|) r 

The terms ±Cm and C y m in (3.24) are integration constants. These factors are obtained 
by considering the u — > oo limit of (2.19), (2.28) and (3.24). Proceeding as in [46], one 
obtains C = — 1 and C y = 0. We have explicitly checked that this procedure yields the same 
integration constant for all possible combinations of the boundary parameter values, namely 
all four regions given in (2.29). 



4 Boundary SSG model and the continuum limit 

To make the paper relatively self-contained, we first give a brief review of the boundary SSG 
model (mainly reproduced from [41]). This is followed by the treatment of the continuum 
limit of the NLIEs given in (3.24). 
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4.1 The boundary SSG model 



The Euclidean-space action of the boundary SSG model is given by 

/oo /»0 /»oo 

dy dx C + / dyCb, (4.1) 
■oo J —oo J —oo 

where the bulk Lagrangian density is given by 

C = 2d z (pd 2 (p - 2^d z tp + 1i\)d- z i\) - 4 cos tp - 4ipip cos ^ , (4.2) 

In (4.2), if) and ■?/> are the two components of a Majorana Fermion field, and z = x + iy, 
z = x — iy. The boundary Lagrangian at x = is given by 

C b = + ia9 y a - 2p((p)a(ip -$) + B((p) , (4.3) 

where a is a Hermitian Fermionic boundary degree of freedom. The functions p((p) and 
B((p), which are potentials that are functions of the scalar field ip, are determined from the 
requirement of boundary integrability and supersymmetry. They are found to be [41] 



1 yF 1 D v sin — 

B(<p) = 2v cos -(<p-<p ), p(<p) = — sin - (tp - D) where tan — = - — — 2 _ - (4.4) 

where F = y/v 2 — 8v cos ^ + 16 

The parameters v and ip are arbitrary and real. As mentioned in [46], in the limit that 
the boundary mass parameters tend to infinity, one arrives at the SSG model with Dirichlet 
boundary conditions. This model corresponds to the open spin-1 XXZ chain with diagonal 
boundary terms [57]. This case was considered in detail in [46]. 



4.2 The continuum limit 

Next, we proceed by analogy with the boundary SG case with general integrable boundary 
conditions [14] and the boundary SSG case with Dirichlet boundary conditions [46]. There, 
in the continuum limit, which consists of taking A — > oo, iV — > oo and lattice spacing A — > 0, 
such that the interval length L and the soliton mass m are given by 

2 -a 

L = NA, m = —e~—, (4.5) 

respectively, the NLIEs of the inhomogeneous open XXZ spin-1/2 and spin-1 chains are 
shown to describe the NLIEs of the boundary SG and SSG models respectively. It is therefore 
natural to conjecture that a set of NLIEs describing a boundary SSG model that corresponds 
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to NLIEs of the inhomogeneous open XXZ spin-1 chain with general integrable boundary 
terms can also be derived here. In the continuum limit, the term — i2iVtan _1 
becomes —i2mL sinh 9 after defining the renormalized rapidity 9 as 



sinh 2^ 
cosh^A 



e = ™. (4.6) 



Thus, (3.24) becomes 

/OO /»oo 
d0' G(0 -9' - ie) ln(l + b(0' + is)) - / dff G(9 -9' + ie) ln(l + b(0' - is)) 
-oo J —oo 

/OO 
d0' G 2 (0 - 9' + ie) ln(l + y(0' - is)) + i2mL sinh 9 + i P Mn/ (0) - ire , 
-oo 

/oo poo 
d9' G(9 -9' - ie) ln(l + b(0' + ie)) + / dff G(9 -9' + ie) ln(l + 6(0' - ie)) 
■oo J — oo 

/oo 
cZ0' G 2 (0' - ^ + ie) ln(l + y(0' + ie)) - i2mL sinh 9 - i P bdry (9) + in , 
-oo 

/oo />oo 
dff G 2 (0 -9' + ie) ln(l + b(0' - ie)) + / dff G 2 {9' -9 + ie) ln(l + b(0' + ie)) 
-oo J —oo 

+ iPytf). (4.7) 



where following definitions have been used, 



e=~, b(0) = &(-), m =!/(-) (4-8) 

/X 7T 7T 



Moreover, G(0) and G 2 (0) are given by 



G(0) = 

7T 7T 



2tt 2 



/oo 
rfA; e" ifc ^ /7r G(ifc) , (4.9) 
-oo 



g 2 (0) = ^g 2 (^: 

7T 7T 



2tt 2 



„* e "**' ft (') = S5E?' (410) 



where G(k) and G 2 {k) are as given in (3.10) and (3.11) respectively. From (3.24), we also 
note that Pbdry{9) and P y (9) are given by 

Pbdry(9) = Pbdry(-) 
71 



JJ_ 

4tt 2 



oo 



dff I dk e - ik ^' /n R(k), (4.11) 
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and 



Py(0) = PyA 
71 



[9 Q 

= 4tt / d6' G 2 (-6') = -2i In tanh- - 2?r , (4.12) 

J -oo 2 

respectively. Note that in (4.11), R(k) is given by (3.25). The boundary terms Pbd ry {0) 
and Py{0) A are essential in the investigations of IR limit (when computing the boundary S 
matrix) and the UV limit (for the computation of effective central charge which is treated 
in section 6). Thus, in particular, PbdryiP) is among the important results of this paper. 



5 Boundary and Casimir energies 

In this section, we compute the boundary correction (order 1) and the Casimir correction 
(order 1/L) to the energy. Following the prescription given in [52] (see also [46]), according 
to which the energy for the inhomogeneous case (A ^ 0) is given by 



E = -i-\ 4-\nA 2 (u) 
A da Iy 1 



where g is given by 



Using (2.28) and the following, 



l/J, 



dk 



[\nh(u)}' = I ^ Lh'(k) e 
2n 



-iku 



ueC 2 , 



(5.1) reduces to 



of dk _,j,a / j±k ffc \ 
E = -4- / — e tkA (e 2 - e a ) 
A ./ 2tt v ; 



Ly'(k)+Lf'(k) + 



27l^J 2 (k) 

ik 



where ip 2 {k) is as defined in section 3. We have also used the following result, 

1 



Lh'(k) 



-ik) 



Lh"(k) 



(5.1) 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



By evaluating the Fourier transform of (2.14) and using (3.21) together with (5.5), we can 
write the energy E as, 



E = El + Ei + Ei/l 



(5.6) 



4 The function P y (0) is the same as for the Dirichlet boundary conditions case found in [46]. 
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where the E L , E-y and Ey/ L represent the bulk vacuum energy, boundary vacuum energy and 
Casimir energy respectively: 



2Nq f°° 9 A , sinh ^ cosh ((^ - %)k) 

E L = — 1 / dke~ 2lAh 2 - ^ V-L . 5.7 

iA /_ m sinh^ 



2 



i A i-oo V 2 / I ch ^ sinh f 



+ + + ( ,S) 



E -JL [°° ^ e -Afc L_ 

1/L " 2Aj_ 00 2n 6 coshf 



LB'(k) + CB'{k) - e- k "LB'(k) - e k »£B'(k)j . (5.9) 



In (5.8), the symbol (+ -H- — ) represents the terms with A + — > A_ , B + — > B_ , s + — > s_. In 
(5.7), adopting the renormalization procedure [29]-[31] of keeping only the finite terms and 
by closing the integral contour in the lower half plane, selecting only the contribution from 
the residue at k = — i^, we note that the integrand is analytic on and within the contour in 
the lower plane, hence giving 

E L = 0, (5.10) 

This is in agreement with known results (see, e.g., [44]). Also see [46]. As found in [46], for 
the case with Dirichlet boundary conditions, next, by using the same contour, we find that 
the boundary vacuum energy (5.8), is given by 

E 1 = m. (5.11) 

where m is given by (4.5). Hence each boundary contributes the vacuum energy m/2, and it 
is independent of the boundary parameters. Finally, the Casimir energy (5.9) simplifies to 

E 1/L = -3m / dd sinh(fl - is) ln(l + b(0 - is)) . (5.12) 

2tt y.oo 

6 UV limit 

We now consider the Casimir energy in the UV limit mL — > 0. A systematic approach 
utilizing dilogarithm trick [47, 53, 54, 55] was used in [46] for the SSG models with Dirichlet 
boundary conditions. The analysis of our set of equations is similar to the one presented 
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there. Thus, we give only the significant differences with respect to [46]. Denoting the 
reference [46] here by I, the Pbdry(oo), (15.5) becomes, 



bdry 



OO 



-R(0) 



= 7T 1 + 



%(s + + s_) _ fi(A + + A_) 
4(f -y) 4(f -/x) 



(6.1) 



The quantity (15.15) becomes, 

2 [| - G(oo)] 



7T (s+ + S_) - /i (A, 



(6.2) 



where G(oo) = 2( ^~ 3 2 / ^ ) ■ In (6.2), u is required to obey the bound \u\ < ^. We further 
remark that the regions of the boundary parameters A± (2.29), should be restricted further 
as follows in order for this bound to be satisfied: 



II, IV 
III 



g<(A + + A_)<| 



-g<(^ + + ^-)<t 
-|< ( A + + A_)<-| 



(6.3) 



The Casimir energy (15.17) is thus given by 



7T 



LE 1/L (0) = -- 



( ir ( s+ + s _)- fi (A + + A.)f 



2 7T(7T-2/x) 

The above result can be rewritten in terms of the effective central charge c e //(0), 



E 1/L (0) 



24L 



'eff 



(0), 



Hence, the result (6.4) implies that the effective central charge (15.18) becomes 

3 3 



"eff 



(0) 



[k(s + + s_)-^A + + A_)Y , 



2 7r(7r-2/i) 

which after imposing constraints (2.10), with k — 1, becomes 



(6.4) 



(6.5) 



(6.6) 



c e //(0) = - 



2 7r(7T - 2/i) 



[vr(s + + s_) =F/x|c_ - c+| 



(6.7) 



We note that the terms in square bracket in (6.6) and (6.7) also appear in the SG case 
with general integrable boundary conditions (refer to equations (2.33) and (2.34) in [14]). We 
also note the similarity with the results obtained in [46], where the boundary SSG model with 
Dirichlet boundary conditions reduces to a system of one free Boson and one free Majorana 
fermion, the central charge of which is given by [58], c = cb + cf — 1 + \ — §■ 
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7 Discussion 



Starting from the T-Q equations for the spin-1 XXZ chain, we derived the NLIEs for inho- 
mogeneous spin-1 XXZ chain with general integrable boundary terms, satisfying a certain 
constraint (2.8). We have restricted our analysis here to the ground state. Further, taking 
the continuum limit, and working in analogy with the usual SG model and the SSG model 
on an interval (with Dirichlet boundary conditions), we also obtain the NLIEs that should 
describe the corresponding SSG model. We further compute the energies that include bulk, 
boundary and Casimir terms. As expected, the bulk contribution vanishes. The boundary 
energy turns out to be independent of the lattice boundary parameters, (A± , B±) (as for 
the case with Dirichlet boundary conditions). We also looked at the UV limit of the Casimir 
energy and obtained an expression for the effective central charge (6.6), that has the right 
c = 3/2 factor as expected. 

One could also further investigate this model. For example, the IR limit can be studied. 
The result obtained for Pbdry{0) is useful for such an analysis. In addition, one could also 
investigate the excited states of the model. Last but not least, it should also be possible 
to carry out such analysis for open spin-1 XXZ chain with non-diagonal boundary terms, 
the solutions to which are given in [56]. In contrast to the solution used in this paper, 
the solutions given in [56] are not restricted by any type of constraints among the lattice 
boundary parameters. It would be interesting to work out the NLIEs for these cases as well 
and analyze their UV and IR limits. We hope to address some of these issues in future. 
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